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Abstract 

We consider nanowires in the field effect transistor device configuration. Modeling each nanowire as a one dimen¬ 
sional lattice with random site potentials, we study the heat exchanges between the nano wire electrons and the sub¬ 
strate phonons, when electron transport is due to phonon-assisted hops between localized states. Shifting the nanowire 
conduction band with a metallic gate induces different behaviors. When the Fermi potential is located near the band 
center, a bias voltage gives rise to small local heat exchanges which fluctuate randomly along the nano wire. When 
it is located near one of the band edges, the bias voltage yields heat currents which flow mainly from the substrate 
towards the nanowire near one boundary of the nanowire, and in the opposite direction near the other boundary. This 
opens interesting perspectives for heat management at submicron scales: Arrays of parallel gated nanowires could be 
used for a field control of phonon emission/absorption. 
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In the elastic coherent regime, thermoelectric trans¬ 
port through a nanosystem involves electric and heat 
currents between at least two electrodes, the source and 
the drain. One has a two-terminal setup where the elec¬ 
trons can be thermalized at different temperatures and 
chemical potentials in each electrodes. In the inelastic 
activated regime, a third terminal must be introduced, 
which contains the (quasi)-particles with which the 
electrons in the nanosystem interact. One often refers to 
this third terminal as the environment. For photovoltaic 
effects, it contains photons. For Mott variable range 
hopping (VRH) regime, it contains phonons. Such 
multi-terminal setups are characterized by two Fermi- 
Dirac distributions describing the thermalized source 
and drain, and a Bose-Einstein distribution describing 
the thermalized baths of photons or phonons. Markus 
Biittiker was known for having set the foundation of 
the theory of multi-terminal mesoscopic transport 1H 
when he was at IBM Yorktown Heights. More recently, 
he studied multiterminal thermoelectric transport. For 
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instance, he considered the case where a third elec¬ 
tron reservoir is capacitively coupled with a nanosys¬ 
tem contacting a source and a drain. In Ref. 0, it was 
shown that an electrical current through the nanosystem 
can cool the third terminal or that heating the third ter¬ 
minal can give rise to an electrical current through the 
nanosystem. This last effect was recently measured in 
Saclay na and Wurzburg ED In Ref. 0 the consid¬ 
ered setup involves three Fermi-Dirac distributions and 
three currents: the electric and heat currents between 
the source and the drain, and the heat current between 
the electrons inside the nanosystem and those inside the 
third terminal. As emphasized |5] (6) by Yoseph Imry, 
hopping transport is also a genuine three-terminal trans¬ 
port involving two electron baths and a phonon bath. 
This is why we feel appropriate to publish this study of 
multi-terminal activated transport in a volume in mem¬ 
ory of Markus Biittiker. 

We consider the array of one dimensional (Id) 
metal-oxide-semiconductor-field-effect-transistors 
(MOSFETs) shown in Fig. [I] The array is made 
of many thin doped semiconductor nano wires (NWs) 
arranged in parallel. Activated transport in each of these 
NWs corresponds to the multi-terminal setup sketched 
in Fig. [2] The source and the drain are the two electron 
reservoirs, whereas the substrate, divided between its 
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source and drain sides, plays the role of a phonon 
reservoir. The electrons are thermalized in the source 
and the drain with Fermi-Dirac distribution of temper¬ 
ature Ts,d and chemical potential jis,d • The phonons 
are thermalized in the substrate with a Bose-Einstein 
distribution of temperature T p h. If the electron states 
are localized inside the NWs, but also coupled with the 
phonons of the substrate, electron transport through 
the NWs becomes activated when their lengths exceed 
the Mott hopping length L M and is governed by Mott’s 
variable range hopping (VRH) mechanism. Hereafter, 
we set jdD = jj = 0 and /i s = b/j = eSV where e is the 
electron charge and SV the bias voltage. We assume a 
uniform temperature Ts = T D = T p h. The bias induces 
an electron current I e from the source to the drain. In 
an activated regime, I e will be also due to the presence 
of local heat currents I Q (x) between the nanowire and 
the substrate. If the nanowire density of states v(E) and 
localization length £(E) do not depend on the energy 
E in the vicinity of yu, I Q (x) fluctuates randomly as a 
function of the coordinate v along the NW, a phonon 
of the substrate being sometimes absorbed, sometimes 
emitted. As a consequence, the ensemble average value 
of the heat exchange vanishes, < I Q (x) >= 0. Interest¬ 
ing effects occur instead if v(E) and £(E) are strongly 
energy dependent. Particle-hole symmetry is broken: 
This gives rise to large thermoelectric effects mm and 
the local heat currents along the NWs become finite, 
< I Q (x) 0, the phonons being mainly absorbed or 
emitted near the NW boundaries. 

We study such a case using for the density of states per 
unit length v(E) and for the localization length %(E) 
analytical results describing a Id lattice with random 
site potentials in the weak disorder limit. Shifting the 
NW conduction band with the gate voltage V g such that 
the electrons are injected from the source around its 
lower edge, the phonons are mainly absorbed near the 
source and emitted near the drain. This corresponds to 
the sketch given in Fig. [2] Inducing with Sjn an electric 
current I e , one obtains a heat current /£ wrce > 0 cooling 
the source side of the substrate, while /£ . < 0 heats 
the drain side. The cooling/heating effects are reversed 
if the electrons are injected around the upper band 
edge and vanish if the injection is made around the 
band center. Phonon absorption/emission can thus be 
controlled with bias and gate voltages in Id MOSFETs. 

1. VRH transport with energy independent localiza¬ 
tion length £ and density of states v 

Phonon assisted hops between localized states and 
the associated mechanism of VRH transport have been 



Figure 1: Thin doped semiconductor NWs (red) are deposited in par¬ 
allel on an insulating substrate (light grey). Two metallic electrodes 
(yellow) provide a source and a drain. An electron current flows 
through the NWs when a bias voltage is applied, while the NW con¬ 
duction bands are shifted by a voltage applied on a back gate (darker 
grey) put below the substrate. The substrate is cut in two parts, the 
middle of the NWs of length L being suspended over a small length 
L c . This defines the source and drain sides of the substrate. 
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Figure 2: Multi-terminal setup corresponding to phonon-assisted 
transport in a single nanowire. It involves two electron reservoirs 
(source and drain) and a phonon bath (provided for a Id MOSFET 
by the insulating oxide substrate on which the NW is deposited). The 
substrate is divided by a small cut between its source and drain sides. 
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mainly studied when the energy dependence of the den¬ 
sity of the localized states per unit length v(E) and of 
their localization length f(E) can be neglected around 
p. Under this assumption, the electron transfer from 
one localized state to another one separated by a dis¬ 
tance x and an energy SE oc 1 /(yx) results from the 
competition between the elastic tunneling probability 
(oc exp{-2v/£}) to do a hop of length x in space and 
the Boltzmann probability (oc exp {-SE/k B T}) to do a 
hop of SE in energy. This competition gives rise to an 
optimal hopping length L M , the Mott’s hopping length, 
and an optimal hopping energy A M . In one dimension, 
one has 
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Electronic transport is achieved via several hops of mag¬ 
nitude ~ L m (with f < L m < L) in space or A M in en¬ 
ergy, and the conductance G can be expressed in terms 
either of L M or A M : 


G oc exp 




(3) 


Since L M is a decreasing function of the temperature, 
the VRH regime takes place above the activation tem¬ 
perature 
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at which L M ^ L and below the Mott temperature 
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at which L M - Below T x , L M exceeds the system 
size and electron transport becomes elastic and coherent 
(see Ref. (3). Above T M , E M < f and transport is sim¬ 
ply activated between nearest neighbor localized states. 
Actually, the crossover from VRH to simply activated 
transport takes place at temperatures lower than T M in 
one dimension. This is due to the presence of highly 
resistive regions in energy-position space, where elec¬ 
trons cannot find available states at distances ~ A M , L M - 
These regions behave as ’’breaks” that electrons are con¬ 
straint to cross by thermal activation, resulting in a sim¬ 
ply activated temperature dependence of the overall re¬ 
sistance idol Ull. 

A more microscopic approach for deriving the above 
expressions consists in replacing the transport problem 
by a random-resistor network one llT2) in which the 


hopping between two localized states i and j is effec¬ 
tively described by a resistor p,y. As proposed by Ambe- 
gaokar, Halperin and Langer, one can approximate 03 
the resistance of the network by the lowest resistance 
p c such that the resistors pp < p c form a percolating 
network traversing the sample. In one dimension, an 
approximation for p c was derived in Ref. ld4l 



which predicts a T~ l/2 behavior independently of the 
length L when T —> 0, and a simple activated ( T~ l ) 
behavior for wires longer than £/2 y/T M /T exp {T M /T}. 


2. VRH transport in realistic thin nanowires 

2.7. Energy dependence of the localization length f(E) 
and density of states v(E) 

All the previous results were based on the assumption 
that f(E) and v(E) are independent of the energy. Such 
an assumption becomes questionable in one dimension, 
notably as one approaches the band edges. Let us take 
the widely used Anderson model: A one dimensional 
lattice of N = L/a sites with lattice spacing a = 1, 
nearest neighbor hopping t, and random site potentials 
6i uniformly distributed in the interval [-W/ 2, W/ 2]. Its 
energy levels E t are distributed in the energy interval 
[E~,E+l, where Ey = ±(2 1 + W/2) defines the edges 
of the NW conduction band when L —> oo. In the bulk 
of the band (i.e. for energies \E\ < 1.50, the density 
v(E) of the Anderson model for small values of W and 
L —> oo can be described by the formula (valid without 
disorder) 

v h (E) = - 1 . (7) 

2?rf V 1 - (£'/2/) 2 


As one approaches the edges E±, the disorder effects 
cannot be neglected and v(E) is given by the analytical 
formula obtained EG) by Derrida and Gardner around 
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where X = (\E\ - 2t)t l/3 (l2/W 2 ) 2/3 and 
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Similarly, the localization length f(E) of the Anderson 
model can be described by analytical expressions valid 
for weak disorder. Inside the bulk of the band, one gets 

^ = $K- £2 )’ ao) 
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Figure 3: Mott hopping length Lm = yJjj/(2vT) as a function of E 
for different temperatures (from top to bottom, kgT = 0.0It, 0.11, t 
and 1(h). v(E) and g(E) are taken from the Id Anderson model. NW 
length L = 1500 and W = t. 



whereas 


&(£) = 2 


/ 12f 2 \ 1/3 

Vw 5 "/ 


j-im 


on 


as E approaches the band edges ±2 l Figures of the 
functions v(E) and %(E) can be found in Ref. m, where 
the weak disorder expansions are shown to fit the nu¬ 
merically calculated values of v(E) and f(E) for values 
of the disorder parameter W as large as t. 


Figure 4: Mott hopping energy A m - kg V TTm as a function of E. 
Same parameters and color code as in Fig. [3] 


2.2. Energy dependence of the VRH scales 

Using the functions v(E) and g(E) of the Id Anderson 
model and Eq. 0,Eq. ([2]) and Eq. ([6]), we show how 
the Mott hopping length L M (E) (Fig. [3]) and the hop¬ 
ping energy A M (E) (Fig. 0 vary with E. The energy 
dependence of the resistance p c which characterizes the 
best percolation network describing the random resistor 
network is shown in Fig. [5] 

Strictly speaking, these functions L M (E ), A M (E) and 
p c (E) do not characterize hopping transport in a sin¬ 
gle Anderson model, but in an ensemble of Id mod¬ 
els where v and f do not depend on E , but rather take 
uniform values which are given by those of the Ander¬ 
son model at an energy E. Nevertheless, one can easily 
guess from the figures of L M (E ), A M (E) and p c (E ) how 
an electron will cross the nanowire. If it is injected from 
the source electrode at an energy around the center of 
the conduction band, it should exchange much less en¬ 
ergy with the substrate (see Fig. [4]) and will find much 
smaller resistances p c (see Fig. [5]) than if it is injected at 
an energy near the NW band edges. Let us take an elec¬ 
tron injected near the lower band edge. It is likely than 
the percolation path that it will follow will consist first 



Figure 5: Logarithm of the resistance lnp c given by Eq. j6| as a func¬ 
tion of E. Same parameters and color code as in Fig. [3] 
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in absorbing many phonons on the NW side close to the 
source, to reach higher energies where the resistance p c 
is weaker (see Fig. |5j. Then, it will continue around 
an optimum energy where it will jump from one local¬ 
ized state to another by emitting or absorbing phonons 
at random. When it will reach the NW drain side, it will 
mostly emit phonons, in order to reach an energy of the 
order of the chemical potential of the drain electrode. 

3. Numerical study of the random resistor network 


use the exact expression (12] ) and we will take for y t j 

Jij = y e pA(gi,€j,rij), (15) 

instead of Eq. ( [T3] ), where 
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To proceed further, we numerically solve the Miller- 
Abrahams resistor network [fT2l which describes hop¬ 
ping transport. The nodes of the network are given by 
the NW localized states. Each pair of nodes i, j is con¬ 
nected by an effective resistor, which depends on the 
transition rates Uy, T induced by local electron-phonon 
interactions. For a pair of localized states i and j of en¬ 
ergies Ei and Ej , Fermi golden rule 0 gives: 

r ij = Jijfi (1 - fj) [Nij + 0{Ei - Ej)] , (12) 

where is the average occupation number of state i and 

= [exp{|£ 7 - - Ei\/k B T} - l] -1 is the phonon Bose 
distribution at energy \Ej - E t \. The presence of the 
Heaviside function accounts for the difference between 
phonon absorption and emission H3l . y t j is the hopping 
probability i —> j due to the absorption/emission of one 
phonon when i is occupied and j is empty. Neglecting 
the energy dependence of £ and v, one obtains in the 
limit Xij » h 

Jij - yep exp(-2 Xij/f ). (13) 

Here x tj = \x t - xj\ is the distance between the states, 
whereas y ep , containing the electron-phonon interaction 
matrix element, depends on the electron-phonon cou¬ 
pling strength and the phonon density of states. Since it 
is weakly dependent on E t , Ej and x t j compared to the 
exponential factors, it is assumed to be constant. From 
hereon we will take y ep = t/h independent of the posi¬ 
tion along the nanowire Under the widely used ap¬ 
proximation [|T3l 1T61 1T71 fT8l \Eij\ » k B T, Eq. (12] ) re¬ 
duces to: 

r ~ y e -2xij/Z e -(\Ei- B \+\Ej- B \+\Ei-Ej\)/2k B T _ ^ 

However, we have seen that in one dimension, the en¬ 
ergy dependence of £ and v cannot be neglected. Here¬ 
after, instead of using the approximation ( fl4| ), we will 


3 This is justified if the length of the cut L c is very small (see Fig. 
This limit will be assumed to hold in all the numerical simulations 
that follow. 


Eq. ( p~5] ) takes into account the energy dependence of 
v(E) and £(E) and is derived in Ref. Q. It reduces to 
Eq.© when & —> gj = £ and (to leading order) when 
& = €»$■ 

The transition rates between each state i and the con¬ 
tacts a (a = L for the source or R for the drain) are 
assumed do be dominated by elastic tunneling contribu¬ 
tions (see Refs. Eii) and read: 

r,„ = y,a f, [1 -fa(E,)\, (16) 

where 

r ia - y e exp(-2 Xia/fr) . (17) 

In the above equations f a (E) = [expKF 1 - /i a )/k B T] + 
1] _1 is the contact cc s Fermi-Dirac distribution, x ia de¬ 
notes the distance of the state i from a , and y e is a rate 
quantifying the coupling between the localized states 
and the contacts. In the following numerical simula¬ 
tions, we will assume y e = t/h. 

Then, the net electric currents flowing between each 
pair of localized states and between states and contacts 
read 

lij = e (Tij - Tji), (18a) 

ha = e (T ia - T ai ), a- L,R (18b) 

e < 0 being the electron charge. Imposing current 
conservation through all the L nodes i of the network 
Yjj hj + Yja ha = 0, one obtains L coupled linear equa¬ 
tions. Solving numerically this set of equations in the 
limit where the bias voltage —> 0 and where the tem¬ 
perature is T everywhere gives the unknown occupation 
numbers of the localized states i and hence the linear 
response solution of the random resistor network prob¬ 
lem (see for more details Refs. 013 ). This allows us to 
obtain all the electrical currents hj and I ia - 

The set of energies E t , of positions x t and of local¬ 
ization lengths hi are required as input parameters of the 
random resistor network problem. Hereafter, we use a 
simplified model as it is conventional in numerical sim¬ 
ulations of VRH transport (see (HJIIfk 0): The E t are 
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uncorrelated variables taken with a distribution corre¬ 
sponding to the density v(E) of the Id Anderson model 
[Eqs. 0 and @], while the corresponding localization 
lengths £( E t ) are taken using Eqs. (TO] ) and (JTTJ). The 
positions x t are taken at random (with a uniform distri¬ 
bution) between 0 and L along a chain of length L = Na , 
N being the number of sites with spacing a = 1. 

4. Heat exchange between substrate phonons and 

NW electrons 

4.1. Local heat currents 

Let us consider a pair of localized states i and j, of re¬ 
spective energies E t and Ej. The heat current absorbed 
from (or released to) the substrate phonon bath by an 
electron in the transition i —> j reads if. = {Ej - Efj if, 
if = Yij - T jt being the hopping particle current be¬ 
tween i and /.I71 The local heat current associated to a 
given localized state i is given by summing over all the 
possible hops with the other NW states j: 

= = (19) 

j j 

We take the convention that if is positive (negative) 
when it enters (leaves) the NW at site i. The probability 
distribution of if is shown in Fig. ^ for states i located 
near three different positions xq along the NW. The gate 
potential (V g = 2.25 1) has been chosen such that the 
chemical potential /i probes the lower band edge of the 
NW conduction band. Near the source, the distribution 
is highly asymmetric: The electrons injected in an en¬ 
ergy interval of width ksT = 0 At around the lower band 
edge need to absorb phonons for reaching higher ener¬ 
gies E where the paths of the random resistor network 
are less resistive. Far from the NW boundaries, the dis¬ 
tribution becomes symmetric: The electrons have found 
an optimum value of E around which they stay, emit¬ 
ting or absorbing phonons at random. Near the drain, 
the distribution becomes again asymmetric, as the elec¬ 
trons emit more phonons to decrease their energies to 
reach the chemical potential // of the drain. 

In Fig. [7] we show 2d histograms of the local heat 
currents as a function of the position x t inside the 
NW. The data have been calculated for a temperature 
ksT = 0.5 1 and four different values of V g , correspond¬ 
ing to electron injection at the band center (V g = 0), 
below the band center (V g = t) and around the lower 
(Vg = 2.250 and upper (Vg = -2.250 band edges of the 
NW conduction band, respectively. At the band cen¬ 
ter, the fluctuations of the local heat currents are sym¬ 
metric around a zero average. They are larger near the 



Figure 6: Probability distribution of if = Ej if f° r 3 positions xo 
along a nanowire [xo = 10 (left), 750 (middle) and 1491 (right)] of 
length L - 1500. Statistics over 2000 NWs, obtained by counting 
the number of ’’points” in [xo - dx/ 2, xo + dx/2], with dx = 20. Pa¬ 
rameters: V g = 2.25 1, ksT = 0 At, Anderson model with W = t, 
6n = 10 _5 r. 



Figure 7: 2d histograms giving the distribution of local heat currents 
if = E j if as a function of the position x ; - along the NW, calculated 
for a gate potential of value V g = 0 (a), V g = t (b), V g = 2.25 1 (c), 
and V g = -2.25 1 (d). Parameters: L = 1500, W = t, ksT = 0.5f, 
6/u = 10" 5 f, statistics over 500 NWs. 
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Figure 8: Mean value < if > of the local heat currents if as a func¬ 
tion of the position x; along the nanowire, for different values of V g 
(V g = 0 (black line), t (green line), 1.5 1 (blue line), 21 (red dashed 
line), and 2.25t (thick brown line)). Parameters: L = 1500 , W = t, 
JzbT = 0.1 1, 6/u = I0~ 5 t, statistics over 2000 wires. Note that for each 
position x/, we count the number of ’’points” in [x/ - dx/2, x; + dx/2], 
with dx = 20 here. 


NW boundaries and remain independent of the coordi¬ 
nate Xi otherwise. Away from the band center, one can 
see that the fluctuations are not symmetric near the NW 
boundaries, though they become symmetric again far 
from the boundaries. When the electrons are injected 
through the NW in the lower energy part of the NW 
band, more phonons of the substrate are absorbed than 
emitted near the source electrode. The effect is reversed 
when the electrons are injected in the higher energy part 
of the band (by taking a negative gate potential V g ): It 
is now near the drain that the phonons are mainly ab¬ 
sorbed. The 2d histograms corresponding to V^ = 2.25* 
and -2.25* are symmetric by inversion with respect to 

if = 0. 

In Fig. [8] one can see how varies the mean value of 
the local heat current if along the nanowire for different 
values of V g . Though the heat current extracted from the 
substrate near the source takes larger local values when 
one probes the band edges, one can see that it is more 
advantageous for cooling the source side of the substrate 
to inject the electrons at an intermediate energy, when 
V g ~ 1.5* (see also Sec. 4.3). 

In Fig. [9j one can see how the standard deviation of 
the local heat currents if varies as a function of the po¬ 
sition Xi along the nanowire, using the same parameters 
and color code than in Fig. [5] The fluctuations of the 
local heat currents are strongly reduced when the elec- 


Figure 9: Standard deviation of the distribution of if as a function of 
the position Xj, obtained using the same parameters than in Fig. [8] 


trons are injected near the band edges, and become in¬ 
dependent of the injection energy if this energy is taken 
around the band center. 


4.2. Total heat current extracted from the source side of 
the substrate 


As sketched in Fig. [T] we assume that the substrate 
has a small cut over which the NW is suspended, and 
we are interested in cooling the substrate on one side of 
the cut, while the other side would be heated. Assuming 
that the contribution of the NW center can be neglected 
(or that the cut is very short), the total heat current ex¬ 
tracted from the source side of the substrate using a NW 
of length L is defined as 

('Ll 2 

Source = < I? > dx, (20) 

Jo 


urce gives the area under the curves shown in 
taken from 0 to L/2. lfource/k B T is shown as 


i.e. If 
Fig- [sj 

a function of V g for k B T = 0.1 1 and 0.5* in Fig. [TO 


One can see that the cooling effect is maximum when 
V g ~ 1.5*. Very approximately, lf ourC e 00 k B T in the 
studied temperature domain and does not depend on L. 
The existence of an asymptotic regime when L —> ex? is 
confirmed in Fig. ED The average local heat currents do 
not depend on L near the NW boundaries and vanish in 
its middle. 
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Figure 10: Total heat current lf ource extracted from the source side 
of the substrate (divided by k B T) as a function of the gate voltage 
V g . Symbols correspond to: L - 1500 and k B T - 0.1 1 (full circles), 
L = 1500 and k B T = 0.5 1 (empty circles), L = 750 and k B T = 0.1 1 
(triangles), L = 3000 and k B T = 0.1 1 (squares) Parameters: W = t, 
5/i - 10 _5 f, statistics over 2000 wires. 



x. 

i 


Figure 11: Mean value < if > of the local heat currents if as a 
function of the position jc/ along the nanowire, for different lengths 
of the nanowire (L = 750 (blue curves), L = 1500 (red curves), and 
L = 3000 (green curves)), and for two values of V g [V g = t (dashed 
lines) and V g = 2 1 (full lines)]. Parameters: W = t, k B T = 0.1 t, 
6fi - 10 _5 f, statistics over 2000 wires. Data for V g = t have been 
shifted upwards for a better visibility. 


4.3. Optimum condition for phonon absorption by the 
first half of a NW 


In Fig. 10 


value for / 


one can see that V g ~ 1.5 1 gives the largest 
urce . This optimal value comes from the 


competition between two effects. On the one side, in¬ 
jecting electrons of energies E near the lower band edge 
favors phonon absorption. On the other side, at these 
energies the resistance p c of the random resistor net¬ 
work is very large (see Fig. [5]), at least near the NW 
boundaries. This necessarily reduces the particle cur¬ 
rent crossing the NW, and hence the heat current ex¬ 
tracted from the source side of the substrate. For a few 
disordered chains, we have studied the map of currents 
larger than some threshold value which are driven by the 
bias voltage. Around the lower band edge (V g = 2.250, 
it does not seem that the electrons are able to hop high 
enough in energy for reaching the band center where p c 
is minimum. Despite very large sample to sample fluc¬ 
tuations, the currents if. connecting states of positions i 
and j located around the middle of the NW have ener¬ 
gies Ei and Ej mainly located at a distance A E « 1.5 1 
below the band center. V g ~ 1.5 1 could correspond to 
the largest energy distance over which the electrons can 
hop for reaching the band center in the middle of the 
NW. The difficulty of reaching energies far from the 
chemical potential p can be understood by looking at 
the expression ( [14] ) for the transition rates. 


4.4. Order of magnitude of heat fluxes using arrays of 
parallel NWs 

In Fig. [TO] the total heat current integrated over the 
first half of a single NW is given in units of t 2 /h. We 
take the parameters values used in Ref. (8): tlk B ~ 
150K gives t « 12.9meV. This implies for lf ource a max¬ 
imum value around 0.4pW when V g = 1.5 1, T = 75K 
and dp = 10 -5 ^ The used bias dpie corresponds to a 
very small voltage ~ 0.13yuV. Assuming that the non 
linear effects remain negligible by increasing the bias 
to a larger value, say dV = 1.3mV corresponding to 
dp = 0 At, one would have a maximum value « 4nW 
for a single NW of a length « 2.2pm. This length corre¬ 
sponds to a lattice of 700 sites with a spacing of 3.2nm. 

Larger cooling power could be obtained by taking 
very large arrays of parallel NWs, as sketched in Fig. [T] 
for 15 NWs only. Let us take a packing density « 20% 
for a 2d NWs array. This corresponds to a NW diameter 
^10 nm and a NW interspacing of « 40 nm. 2 x 10 5 
parallel NWs would make an array of width » 1cm. It 
could be used for extracting « lmW from the source 
side of the substrate, in a 1 cm large and a few pm 
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long area located near the source electrode. This esti¬ 
mated value of « lmW is done assuming at a temper¬ 
ature T « 75K and a bias SV = 1.3mV. It could be in¬ 
creased by increasing T or SV, if VRH transport holds 
at larger temperatures and if the linear response theory 
remains valid at larger bias voltages. Additional results 
concerning the perspectives given by large arrays of Id 
MOSFETs for energy harvesting and hot spot cooling 
can be found in Ref. 0). 

5. Summary 

We have considered phonon assisted transport be¬ 
tween localized states using a model characterized by 
the density of states v(E) and the localization length 
£(E) of the Id Anderson model, going beyond theories 
where the energy dependence of v and f are neglected. 
We have studied the distributions of the heat currents 
characterizing locally the heat exchange between the 
substrate phonons and the NW electrons, as one shifts 
with a gate the NW conduction band. From those local 
heat currents, the structure of the percolative network 
connecting states of different energies E t and locations 
x t can be guessed. When the electrons are injected in 
the middle of the conduction band, the percolative path 
consists mainly of states with energy E t in an energy 
interval of width A M around the band center. Even in 
this case, Fig. [7] (a) shows us that the heat currents have 
larger fluctuations near the NW boundaries that in its 
bulk. This has to be related to the observation made 
in Ref. (6) that the thermopower in the hopping regime 
is governed by the edges of the samples. When elec¬ 
trons are injected around the band edges, the boundary 
effects upon the local heat currents become much larger. 
It becomes possible to have mainly phonon absorption 
at one NW boundary, and phonon emission at the other, 
instead of spreading phonon emission and absorption 
along the NW. In Refs. (71 El, the thermoelectric effects 
of gated disordered NWs have been studied in the hop¬ 
ping regime. In this work, we have focused the study on 
the heat exchanges occurring between the NW electrons 
and the substrate phonons. The obtained results lead us 
to put forward arrays of Id MOSFETs as tools for man¬ 
aging heat at submicron scales. Moreover, such tools 
make possible to control phonon absorption/emission 
by varying the gate and bias voltages. 
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